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Several  variational  principles  are  derived  for  the  initial¬ 
boundary- value  problem  of  fully  coupled  linear  thermoelasticity  for 
an  inhomogeneous,  anisotropic  continuum.  A  consistent  set  of  field 
variables  is  employed  and  a  method  based  on  the  Laplace  transform 
is  used  to  incorporate  the  initial  conditions  explicitly  into  the  form¬ 
ulation.  These  principles  lend  themselves  readily  to  numerical 
solutions  based  on  an  extended  Ritz  method. 
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1 .  Introduction^ 

The  application  of  variational  methods  for  both  the  unified 
development  of  the  theory  and  for  the  approximate  solution  of  fully 
coupled  initial-boundary- value  problems  in  linear  thermoelasticity  is 
not  new.  As  a  starting  point,  the  work  of  Biot  [l]  introduced  a  varia¬ 
tional  principle  in  terms  of  a  pair  of  vector -valued  primary  variables, 
the  displacement  of  a  material  point  and  a  variable  which  he  termed 
the  entropy  displacement.  The  Euler  equations  of  the  principle  are 
the  thermoelastic  equations  of  motion  and  the  corrected  heat  con¬ 
duction  equation.  A  generalization  of  Biot’s  principle  [2]  incorporated 
the  additional  primary  variables  of  the  stress  tensor  and  a  thermal 
dis -equilibrium  force  conjugate  to  the  entropy  displacement.  Euler 
equations  representing  the  linear  thermoelastic  stress-strain  rela¬ 
tionship,  the  Fourier  heat  conduction  law,  and  a  relationship 
between  the  temperature  and  the  thermal  gradient,  as  well  as  the 
aforementioned  field  equations,  are  products  of  that  variational 
principle . 

The  most  general  variational  statement  of  the  coupled 
thermoelastic  problem  was  made  by  Bao-Lian1  Fu  [3]  and  later  by 
Ben-Amoz  [4],  both  of  whom  obtain  essentially  the  same  Euler 
equations  and  natural  boundary  conditions  Both  prescribe  a  bound¬ 
ary  condition  on  the  entropy  displacement  vector  rather  than  on  the 
heat  flux  vector,  and  Ben-Amoz  obtains  a  sixth  Euler  equation, 
representing  a  relationship  between  the  temperature  and  the  thermal 
gradient. 

Impetus  for  further  development  of  variational  principles  for 
the  coupled  thermoelastic  problem  is  suggested  by  recent  work  of 
Gurtin  [5],  [6],  In  these  treatments  of  linear  elastodynamics  and 
transient  heat  conduction  Gurtin,  utilizing  the  operational  methods 
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of  Mikusinski  [7],  explicitly  introduces  the  initial  conditions  appro¬ 
priate  to  the  problem  into  the  field  equations  and  governing  functionals, 
and  derives  alternate  characterizations  of  the  problems  The  follow¬ 
ing  work  represents  an  extension  of  these  concepts  to  the  field  theory 
of  linear  coupled  thermoelasticity,  and  an  attempt  to  remove  a  lack 
of  consistency  (present  in  the  existing  variational  formulations)  in 
the  choice  of  the  field  variables  and  boundary  conditions 
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2.  Notation^  Mathematical^  jprelimin  arjies^ 


The  notation  and  format  used  by  Gurtm  [5]  will  largely  be 


employed  here.  Therefore  the  standard  mdicial  system,  with  Greek 
and  Latin  subscripts  ranging  over  the  integral  values  (1,2)  and 
(1,2,3)  respectively,  is  used  throughout.  The  reference  frame  is 
Cartesian,  repeated  subscripts  imply  summation,  and  Kronecker's 
delta  is  denoted  by  5...  Parenthetical  superscripts,  as  well  as 


superposed  dots,  indicate  the  order  of  time  differentiation,  and  sub¬ 
scripts  preceded  by  a  comma  denote  space  differentiation  with 
respect  to  the  Cartesian  coordinates.  Parentheses  about  a  pair  of 
free  subscripts  will  signify  the  symmetric  part  of  the  tensor  with 
respect  to  those  subscripts. 


A  region  V  shall  denote  the  closure  of  an  open,  bounded, 


connected  set  contained  in  the  three-dimensional  Euclidean  space  E 
The  boundary  of  V  we  denote  by  S,  being  the  union  of  a  finite  number 
of  non-intersecting  closed  regular  surfaces.  The  term  regular 


surface  is  used  in  the  sense  of  Kellogg  [8]  The  interior  of  V  is  V 
and  n  is  the  outward  unit  normal  vector  to  S.  S  and  J  aie  a 


ai  e  a 


a 


a 


dual  system  of  complementary  regular  subsets  of  S,  where  and 

denote  the  closures  of  and  respectively. 


The  domain  of  definition  for  all  functions  of  the  position 


vector  x  and  the  time  t  is  the  set  which  is  the  Cartesian  product 
of  the  region  of  space  V  and  the  interval  of  time  [0,°°)  denoted  by 
V  x  [0,®). 


The  values  of  a  function  f(x,  t)  and  its  derivatives  are 
defined  on  the  boundary  of  the  domain  of  definition  V  x  [0,°°)  by 


(x  ,  t  )  =  lim 


f  n) 

lim  f'..  (x,  t)  , 

1J.  .  .  K 

(x,  t)  -►  (x  ,  t  ) 

~~o  o 


(2.  1) 
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where  x  eV,  t  =  0  or  x  eS,  t  e[0  °°)  We  say  that  the  function  f 
~o  o  ~o_T  o 

M ,  N  —  f 

is  in  the  function  class  C  if  and  only  if  it  is  defined  on  V  x  [0,  ®), 

and  all  of  the  functions 


exist  and  are  continuous  on  V  x  [  0,  °°) . 

A  point  xeS,  or  a  point  (x,  t)eSx[  0,  ®),  such  that  n  is  con¬ 
tinuous  at  x,  will  be  termed  a  re^guj^r^jgodnt^  A  function  f  will  be 
said  to  be  piecewise  regular  on  S  x  [  0,  »)  if  and  only  if  f  is  piece¬ 


a 


wise  continuous  on  S  x  [  0,  <»),  and  every  regular  point  of  that  region 
is  a  point  of  continuity  of  f.  Two  piecewise  regular  functions  defined 
on  S  x  [0,<»)  are  said  to  be  equal  if  and  only  if  they  are  equal  at 
every  regular  point  (x,  t)  e  x[0,<»). 

Let  f  and  g  be  functions  of  space  and  time  defined  on 
V  x  [0,  <»)  such  that  both  are  continuous  on  [0,  ®)  for  each  x  e  V 
The  of  ^  and  6  given  by 

[f*g]  (X,  t)  =  J^ffx,  t-T)  g(x,  T)dT,  (x,  t)e  Vx[  0,  °°),  (2,3) 


having  the  well-known  properties 


(a)  g*f  =  f*g  (2.4a) 

(b)  g*(f*h)  =  (g*f)*h  =  g*f*h  (2  4b) 

(c)  g*(f+h)  =  g*f+g*h  (2.  4c) 

A  corollary  of  Titchmarsh's  theorem  [9]  will  prove  useful  in  later 
work  and  is  introduced  in  the  form 
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f*g  =  0  implies  either  f=0  or  g  =  0.  (2.  5) 

The  term  functional  will  be  used  to  identify  a  real-valued 
function  whose  domain  is  a  subset  of  a  linear  space.  If  L  is  a 
linear  space,  K  a  subset  of  L,  and  Q^{'}  a  functional  defined  on 
K ,  then  for 

R,  R  e  L, 

R  +  X  R  e  K  for  every  real  X,  (2.  6) 

formally  define  the  notation 


6  OtR}  =  4-  n  [R  +  XR)  I  X  =  0.  (2.  7) 

R 

The  variation  of  Q{  *  }  is  zero  at  R  over  K  and  is  written 

6Q{r]  =  0  over  K,  (2.8) 


if  and  only  if 


6  n{R} 

R 


exists  and  equals  zero  for  every  choice  of  R 


consistent  with  (2.  6). 

Three  lemmas  which  are  analogous  to  the  fundamental  lemma 
of  the  calculus  of  variations,  and  which  have  been  proved  by  Gurtin 
[5],  are  now  stated  Based  on  these  lemmas  we  also  state  a 
corollary  which  will  be  needed  in  the  subsequent  development. 

2.  1  Lemma,  Let  f  be  a  continuous  function  on  V  x  [0,°°), 


and  suppose 


J  [f*g]  (X,  t)  dV  =  0  (0  ^  t  <  oo) 

V 


(2.  9) 
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OO  f  00 

for  every  g  0  C  which,  together  with  its  space  derivatives, 
vanishes  on  S  x  [0,oo).  Then 

f=0onVx[0,°°).  (210) 

L^_2 .  Lemma,  Let  f  be  a  piecewise  regular  function  on 
S2  x  [0,a>)(  and  suppose 


J  [f*g]  (x,  t)dS  =  0  (0  $  t  <  «) 


(2  11) 


00  00  p 

for  every  g  0  C  that  vanishes  on  x  [O,00)  Then 


f  =  0  on  S^xfO,00) 


(2  12) 


2_.^_3 ,  Lemma 

suppose 


Let  f. 

1 


be  continuous  on 


Sl  X  [0,  ®), 


and 


f  [f  *(g..n.)]  (x,t)dS  =  0  (0  ^  t  <  00)  (2.13) 

V  1  1J  J  ~ 

bl 

00  OO 

for  every  g.,  0  C  ’  which,  together  with  all  its  space  derivatives, 
vanishes  on  x  [O,®)  and  which  has  the  property 


(2.  14) 


Then 


f.  =  0  on  S  x  [0,  ® ).  (2.  15) 

1  1 


A  corollary  of  Lemma  2.  3  is; 
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2.4.  Corollary.  Let  f  be  continuous  on  x  [  0 , 00 ) ,  and 

suppose 

J  [f*(g.n. )]  (x,  t)  dS  (0  2  t<  •)  (2  16) 

J1  1  1  ~ 

CD  00 

for  every  g.  e  C  which,  together  with  its  space  derivatives, 
vanishes  on  4  x  [  0,  °°).  Then 

f  =  0  on  J  x  [O.oo).  (2.  17) 
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3.  The  initial-boundary-value  problem  with  mixed  boundary 
conditions . 


The  fundamental  system  of  field  equations  for  linear  coupled 
thermoelasticity,  characterizing  inhomogeneous  and  anisotropic 
solids,  is  now  stated  for  reference.  Let  V  be  the  region  of  space 
occupied  by  the  solid  and  let  V  x  [0,  »)  denote  the  domain  of  defini¬ 
tion  for  all  functions  of  position  and  time. 

Let  u(x,  t),  e..(x,t),  T..(x,t),  F.(x,t),  q.(x,  t),  and  (x,  t), 
i  -  ij  —  ij  —  i  ~~  i  ~~  i  ~~ 

in  this  order,  represent  the  Cartesian  components  of  the  displace¬ 
ment  vector  u(x,  t),  the  infinitesimal  strain  tensor  e^x,  t),  the 
stress  tensor  T(x,  t),  the  body  force  vector  F(x,t),  the  heat  flux 
vector  jq(x,  t),  and  the  thermal  gradient  vector  ${x>  t).  Also  let 
0  (x,  t)  denote  the  temperature  above  a  quiescent  reference  state 
(a  constant  absolute  temperature),  and  let  r|(x,  t)  and  H(x,  t), 
respectively,  be  the  specific  entropy  per  unit  mass  and  the  rate  of 
internal  heat  generation  per  unit  volume  within  the  solid. 

Then  the  strain-displacement  and  thermal  gradient-temperature 
relations  are 


e . 
ij 


=  —  (u.  .  +  u.  .)  = 

2  l.J  J-1 


u 


(i.  j) 


(3.  la) 


and 


i?.  =  9,.  on  V  x[0,°°). 

l  l 


(3,  lb) 


The 


re  spectively. 


T. .  .  +  F.  =  p  ii.,  T..  =  T„., 

iJ.J  i  1  ij  J1 


(3.  2a) 
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and 


q.  +  p  T  f|  =  H  on  V  x  [0,  »), 

1  5  1  O 


(3.  2b) 


where  p(x)  is  the  mass  density  of  the  solid.  Then  introduce 

c...  _(x),  k.  (x),  a  , (x) ,  C  (x),  and  C  (x)  as,  respectively,  the  com- 
ijkl  -  ij-  ij-  e~  T  ~~ 

ponents  of  the  isothermal  elasticity  tensoi  ^(x),  the  thermal  conduc¬ 
tivity  tensor  k(x),  the  thermal  expansion  tensor  oc(x),  and  the 
scalar  specific  heats  for  zero  defoimation  and  for  zero  stress.  The 

three  sets  of  tensor  components,  X..,  ,(x)>  X  (x),  and  8  (x),  will  be 

ijkl  ~  ij  ~  ij  ~ 

termed  the  components  of  the  isothermal  compliance  tensor  X(x), 
the  thermal  resistivity  tensor  XJx),  and  the  thermoelasticity  tensor 
J3(x).  These  quantities  are  related  through 


c.  VI  xw  =  $.  6  , 

ijkl  klmn  1m  jn 


k.  \  .  =  6.., 

lm  mj  ij 


P..  -  c.  a 
ij  ijkl  kl 


C  =  C  -  a.  (3  T/p  on  V, 
e  T  ij  ij  o 


and  satisfy  the  symmetry  conditions 


(3.  3a) 
(3.  3b) 
(3.  3c) 
(3.  3d) 


Cijkl  °jikl  klij 

(3  4a) 

Xijkl  =  Xjikl  “  Xklij’ 

(3,  4b) 

k. .  =  k  . , 

(3  4c) 

ij  J1 

X..  =  X.., 

(3.  4d) 

!J  J1 
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a. .  =  a  , 
iJ  J1 


3  . .  =  3  .  on  V. 
ij  J1 


(3  4e ) 
(3  4f) 


The  constitutive  relations  and  the  equation  of  state  are  then 


and 


T  = 

1J 

°  ijkl  6kl 

q.  = 

-k  9  , 

i 

ij  J 

P  T 

T)  =  P  C  0 

ij  o  lj 


(3  5a) 
(3  5b) 


or,  alternatively, 


e.  =  X.  ,  ,  T  +  a..  0, 
ij  ljkl  kl  ij 


q., 

i  ij  J 


and 


pT  ri=pC  9  +  a. .  T  t.,  on  V  x  [0, »). 

O  T  IJ  O  lj 


(3  6a) 
(3.  6b) 
(3.  6c) 


Associated  with  this  system  of  field  equations  are  the  initial 
conditions 


and 


u.  (x,  0  )  =  d,  (x) , 

l  ~  l  ~ 


u.(x,  0)  =  v  (x), 


0  (x,0)  =9  (x)  on  V, 

~~  O  ~~ 


the  displacement  boundary  conditions 


u  =  u  on  S,  x  [0,  oo), 
ii  1 


(3  7a) 
(3.  7b) 
(3  7c) 

(3.  8a) 
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the  traction  boundary  conditions 


T.  =  T.  n  =  T  on  S_  x  [0,®), 

i  ij  J  i  2 


the  temperature  boundary  condition 


9=9  on  x  [0,  “), 


(3.8b) 


(3  8c) 


and  the  he  at  flux  jp o undar y^jco nd it i o n 

Q  =  q.  n  =  Q  on  J  x  [0,®).  (3  8d) 

ii  L 

Here  d,  v,  and  0  are  the  prescribed  initial  displacements,  initial 
—  o 

velocities,  and  initial  temperature  distribution,  while  u,  0,  and  Q 
are  the  given  surface  displacements,  tractions,  temperature,  and 
normal  heat  flow. 

It  should  be  noted  that  the  boundary  conditions  prescribed 
above  may  be  generalized  to  include  mixed-mixed  conditions  for 
tractions  and  displacements,  as  well  as  the  more  general  ’’radiation’1 
type  boundary  conditions,  such  as  the  elastically  supported  surface 
and  the  thermal  convection  boundary  layer 

Analogous  to  the  work  of  Gurtin  [5],  the  smoothness  require¬ 
ments  and  other  regularity  assumptions  on  the  ascribable  functions 
are  introduced  as  hypotheses  on  the  data' 

(i)  p  >  0  is  continuously  differentiable  on  V; 

(ii)  jc  and  X  are  continuously  differentiable  on  V  and 
satisfy  (3  3a),  (3,4a),  and  (3  4b); 

(iii)  C  >  0  and  C  >  0  are  continuously  differentiable  on 

_e  T 
V  and  satisfy  (3  3d); 


13 


(iv)  k,  a,  X^,  and  are  continuously  differentiable  on  V 
and  satisfy  (3  3b),  (3  3c),  (3  4c),  (3  4d),  (3.  4e), 
and  (3  4f); 

(v)  d  is  continuously  differentiable  on  V; 

(vi)  v  and  9  are  continuous  on  V; 


o 


(vii)  F  and  H  are  continuously  differentiable  on  V  x  [0  °°); 


(viii)  u  and  9  are  continuous  on  x  [0,®)  and  x  [o,®), 
re  spectively ; 


•  •  A  *i  p 

(ix)  T  and  Q  are  piecewise  continuous  on  x  [0,®)  and 
J '  x  [  0  ® ) ,  respectively. 


With  the  specifications  on  the  data  cited  above  the  mixed 


problem  consists  of  finding  a  set  of  functions  [u,  jT,  9,  T\,  d] 


on  V  x  [0,®)  which  satisfies  the  field  equations  (3  1),  (3  2),  (3  5), 
the  initial  conditions  (3  7),  and  the  boundary  conditions  (3.  8). 


The  ultimate  goal  is  to  transform  the  foregoing  statement  of 


the  mixed  initial-boundary -value  problem  into  an  equivalent  varia¬ 
tional  formulation.  To  accomplish  this,  it  is  convenient  to  define  an 
admissijple  state .  and  then,  in  terms  of  admissible  states,  a  solution 


0,  Tl,  3.’  $)•  is  an  ordered  array  of  functions  u,  ^e,  JT,  0,13,  c[,  ^ 
defined  on  V  x  [O,®)  with  the  properties: 


(a)  u  e  C 


i 


1,2 


,  e  .  e  C 


i  i 


03  0 


(3  9) 


on 


V  X  [0,  o=>). 


Addition  of  states  and  multiplication  of  a  state  by  a  scalar  are  defined 
by 
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R  +  R  =  [u+Tt,  e+e,  T+T,  0+S',  n+?f,  £+;§,  £+£] 
and  (3.  10) 

\R  =  [\u,  Xje,  Xt_,  X8,  Xr) ,  Xc[,  \j>]. 


In  this  way  the  set  of  all  admissible  states  is  a  linear  space, 

A  solution  of  the  mixed  problem  is  now  defined  as  an  admis¬ 
sible  state  R  =  [u,  je,  jt,  8,  r\,  c|,  &]  which  satisfies  the  field  equa¬ 
tions  (3.  1),  (3  2),  (3,  5),  the  initial  conditions  (3.  7)  and  the  boundary 
conditions  (3.  8). 
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4.  Aljejr native  f o r ul a _t  1  ons . 

Since  a  major  motive  for  the  recasting  of  the  initial -boundary- 
value  problem  of  coupled  thermoelasticity  is  to  incorporate  the  initial 
conditions  explicitly  into  the  field  equations  and  into  the  functionals 
which  arise  in  the  variational  formulations,  it  is  convenient  to  con¬ 
sider  the  Laplace  transformation  of  equations  (3.  2).  Using  (3  5c) 
and  (3.  7) 


T..  ,  +  F.  +  osd  I  ov,  =  os  u  (4,1a) 

ij ,  J  i  i  i  i 


and 

q  -  pC  0  -  p.  T  d.  +  oT  sr\  =  H,  (4  lb) 

l,  i  e  o  ij  o  i,  j  o 

where  a  superimposed  bar  denotes  the  transformed  function  and  s  is 
the  transformation  parameter.  Solving  (4  1)  for  u.  and  T),  then 
applying  the  inverse  transformation,  yields 


Ou  =  g*T  ..  .  +  f. 

i  ijij  i 


(4  2a) 


and 


pt  n 

o 


h  -  g'*q. 


i,  i 


(4.  2b) 


where  the  functions  g,  g1,  f.,  and  h  are  defined  to  be 


g(t)  =  t,  g'(t)  =  1  (0  5  t<  «),  (4.  3a) 

f.(x,  t)  =  [g*Fj  (x,  t)  +  o(x)  [t  v.(x)  +  d.(x)],  (4,  3b) 


and 
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h(x,  t)  =  [g'*H]  (x,  t)  +  p(x)  C  (x)  8  (x) 


+  T 


0.  (x)  d  (x),  (x,  t)  e  V  x  [ 0 ,  . 

o  ij  ~  i5J 


(4  3c) 


The  functions  f,  and  h  are  completely  described  by  information 
relevant  to  the  mixed  problem 

With  these  results  in  mind  alternative  formulations  of  the 
problem  can  easily  be  made, 

4.1.  Theorem.  Let  ut  £  T. .  e  ,  and  suppose 

—  - -  i  ij 

T..  =  T  „ . .  Then  u  ,  T , .  satisfy  the  equations  of  motion  (3.  2a)  as  well 
J1  1  IJ 

as  the  initial  conditions  (3  7a)  and  (3  7b)  if  and  only  if 


g*T  .  .  +  f  =  pu  on  V  x  [  0,  co).  (4.4) 

ij .  J  1  1 

The  proof  of  this  theorem  has  been  given  by  Gurtin  [  5] 

4^2  Theorem  Let  T)  eC^  \  q  e  and  suppose  the 

equation  of  state  (3.  5c)  holds  for  t  =  0  Then  T]  and  q„  satisfy 
the  energy  equation  (3  2b)  as  well  as  the  initial  conditions  (30  7c)  if 
and  only  if 


h  -  g 1  *  q  .  =  p  T  T|  on  V  x  [  0,  » ) .  (4.5) 

1,1  o 

Proof,  Suppose  (3  2b),  (3.  7c),  and  (3,  5c)  (for  t=0)  hold 
Then  (4.  3a)  implies 

[g'*(H-q.  )]  (x,  t)  =  p(x)  T  J^f)  (x,  T)  dr 

1,1  o  J0 


=  p(x)  T  r)  (x,  t)  -  p(x)  T  T)  (x,  0) 

O  ^  O  ~~ 

=  p(x)  T  r)  (x,t)  -  p ( x)  C  (x)  9  (x) 
o  ~-e~“o~- 

-T  3 .  (x)  d  (X).  (4.6 

o  ij  ~  i,  J  ~ 
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Thus,  and  by  (4.  3c),  equation  (4.  5)  is  met.  Conversely,  suppose 
(4.  5)  holds.  Then,  by  reversing  the  argument,  and  utilizing  the 
definitions  (4.  3a)  and  (4.  3c),  it  is  directly  verified  that  T|,  q  meet 
(3.  2b).  Since  (3.  5c)  (for  t-0),  (4.  3a),  (4.  3c),  (4.  5)  imply  (3.  7c), 
the  proof  is  complete. 

4.3 .  Theorem.  Let  R  =  [u,  _e ,  T^,  0,  Tl,  <£>  j?]  be  an  admis¬ 
sible  state.  Then  R  is  a  solution  to  the  mixed  problem  of  coupled 
thermoelasticity  if  and  only  if  it  meets  the  field  equations  (3.  1), 

(3.  5),  (4.  4),  (4.  5),  and  the  boundary  conditions  (3.  8). 

This  result  is  a  trivial  consequence  of  Theorems  4,  1  and  4.  2. 

As  a  result  of  this  theorem,  an  alternative  characterization 
of  the  solution  to  the  mixed  problem  of  coupled  thermoelasticity  has 
been  developed  such  that  the  initial  conditions  are  explicitly  incor¬ 
porated  into  two  of  the  field  equations. 

By  a  d i Lg_pl a c e mnent  and jt^e mpe r a^ur^ield^corres 
solu^tion_  oi  t he^m i  x  e  d  p  r o bhem  i s  meant  an  ordered  pair  [u,  0]  made 
up  of  a  vector-valued  function  u  and  a  scalar- valued  function  0  such 
that  there  exist  functions  e,  T,  r\,  c[,  $  with  the  property  that 

[u,  je,  JT,  0,  T] ,  c[,  $]  is  a  solution  to  the  mixed  problem. 

2  2  2  0 

4.  4  Theorem.  Let  u„  €  C  and  0  e  C  .  Then  u,  0 
is  a  displacement  and  temperature  field  corresponding  to  a  solution 
of  the  mixed  problem  of  coupled  thermoelasticity  if  and  only  if 


and 


g*(  c  i*|  Ui 

ljkl  k, 


P..0).  •  +  f. 

ij  J  i 


pu 

i 


(4.  7a) 


g'*(k  0.  ),  .  +  h  =  pC  0  +  T  p  u .  on  Vx  [0,®)f 

ij  ij  i  e  o  ij  i,j 

u^  =  u^  on  S  x  [O.oo),  (4.  7c) 


(4.7b) 
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(  c  u  -  P  8  )n.  =  T  on  S  x  [0,°°),  (4.  7d) 

ijkl  ]>1  ij  J  i  2 

9=8  on  ^  x  [0  «),  (4.7e) 

and 

-k  8,  n.  =  Q  on  f  x  [0,«)  (4.  7f) 

ij  J  i  2 


^Proof.  First  sappose  a  and  0  meet  equations  '4.  7}0  Detine 
e^  and  $  through  (3.  1)  and  JT,  jq,  T|  through  (3,  5).  Then,  and  because 
of  (3.4),  the  boundary  conditions  (3.8)  are  met.  Further  (4  7a), 

(4.  7b),  (3.  1),  (3.  4),  (3.  5)  imply  (4  4),  (4  5).  Thus,  and  by  Theorem 
4.  3,  [u,  je,  JT,  0,  p,  jq,  #]  is  a  solution  to  the  mixed  problem, 
Conversely,  (3.  1),  (3.  4),  (3  5),  (3.  8),  (4  4),  (4  5)  imply  (4.  7)  and 
the  proof  is  therefore  complete. 

By  a  s tress _and.  h  to  a  s o lution  of 

the  mix gd  problem  is  meant  an  order  ed  pair  [jT,  jq]  made  up  of  a  sym¬ 
metric  second-order  tensor- valued  function  and  a  ve ctor -valued 
function  jq  such  that  there  exist  functions  ^u,  je,  0,  p,  &  with  the 

property  that  [u,  e,  T,  0,  p,  jq,  ??]  is  a  solution  of  the  mixed  problem. 

~~  2,0  2  0 

4.5.  Theorem,  Let  T  e  C  and  q  e  C  with  T.  = 

™  - -  ij  i  ij 

T  ..  Then  T ,  q  is  a  stress  and  heat  flux  field  corresponding  to  the 
J 

solution  of  the  mixed  problem  of  coupled  thermoelasticity  if  and 
only  if 


g*T 


(im,  m 


(h  -  g'*q 


m, 


) 

m 


v  1  T 

ljkl  kl 


(4  8a) 
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and 


a' .  t 


1 


kl  kl  pC 


(h 


g'*q  )~1  .  =  V  q 

m,  m  J  ,  1  ij  j 


on  V  x  [  0,  «>) , 


(4  8b) 


—  (g*T  .  +  f.)  =  u. 

P  1J,J  i  i 


on  x  [0,=°;,  (4.  8c) 


t  n.  =  T  =  T 
ij  J  i  1 


on  S2  x  [  0,  »),  (4  8d) 


— — -  (h  -  g'*q  )  -  a!  t..  =  9 
PC  B  Hm,  rn  ij  ij 


on  J.  x  [0,<=), 


(4. 8e) 


and 


q.  n  =  Q  =  Q 

l  l 


on  xju  x  [0, °°),  (4.  8f) 


whe  r  e 


T  a.,  (x) 

a.  ,'X)  =  -9 

ij  v~-  P(x)  C  (x) 


X!.,,(x)  =  X  ..(x)  -  a  (x)a'  (x). 

ljkl  ljkl  —  IJ  —  kl  — 


(4  9) 


Proof  First  suppose  jt  and  meet  equations  (4.8), 
Define  u  from  (4.  4),  T|  from  (4  5),  0  from  (3  6c),  £  from 
(3.6a),  and  £  from  (3.  6b).  Then  (4.  8a),  (4.8b),  (4  8c),  (4.  8e), 
(4.9),  (3.6),  (4.4),  (4  5)  imply  (3.  1),  (3,8a),  (3.8c).  Finally  (3,  5) 
hold  by  virtue  of  (3.  3),  (3.  6).  Thus,  and  by  Theorem  4C  3, 

[u,  'H,  jcj,  #]  is  a  solution  of  the  mixed  problem 

Conversely;,  (3.  1),  (3.  3),  (3  5),  (4.  4),  (4,  5),  (3  8)  imply  (4.  8)  and 
the  proof  is  complete. 
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The  fact  that  the  behavior  of  the  stresses  in  linear  elastody- 
namics  can  be  characterized  by  a  single  tensorial  field  equation, 
similar  in  form  to  (4.  8a) ,  is  due  to  Ignaczak  [  1 0] .  Here,  by  (4.  8a) , 
(4.  8b),  this  has  been  extended  to  the  linear  theory  of  dynamic, 
coupled  thermoelasticity. 
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tjierm^ 


The  most  general  variational  principle  to  be  considered  here 


will  be  one  in  which  the  admissible  states  (which  constitute  a  linear 
space)  are  not  required  to  meet  any  of  the  field  equations,  initial 
conditions,  or  boundary  conditions  A  number  of  less  general  prin¬ 
ciples  can  also  be  derived,  depending  on  the  extent  to  which  certain 
of  the  requirements  are  taken  to  be  identically  satisfied  by  the 
admissible  states.  Three  principles  will  be  formulated  here,  coi  re¬ 
sponding  roughly  to  the  elasto- static  variational  principles  of 
Hu-Washizu  [  1 1] ,  [  1 2] ,  Hellinger-R  eis  sne  r  [  1 3] ,  [  14] ,  and  minimum 
potential  energy  [  1 5 J . 

For  consequent  use,  recall  from  equations  (4.  3)  that 


g(t)  -  t,  g'(t)  =  1  ( 0  =  t  <  «), 


(5,  la) 


f.(x,  t)  =  [g*F.J  (x,  t)  +  p(x)  [tv.(x)  +  d  (x)],  (5.  lb) 

1  1  -  -  1  -  1  ~ 


and 


h(x.  t)  =  [g'*H]  (x,  t)  +  p(x)  C  (x)  9  (x) 

—  ~~  e  o  ^ 


(5. 1c) 


+  T  9  (x)  d  .(x) ,  (x,  t)  e  V  x  [0,“). 

o  ij  -  i,J~  ~ 


Also  T  and  T  shall  be  written  for  the  traction  vector  s  with 


components 


T. 


l 


T . .  n  ,,  T.  =  T, .  n 
J  1  J 


(5.  Id) 


and  Q  and  Q  for  the  normal  heat  fluxes  given  by 


Q  =  q.  n.,  Q  =  q,  n.. 
ii  ii 


(5  le) 
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5VL 

Let  R  =  [u, 
functional  C 


ThepRjri,  Let  K  be  the  set  of  all  admissible  states 
je,  JT,  9 >  r|>  J9>  e  K,  and  for  each  t  e[0,»)  define  the 

t(  *  3  on  K  by 


O  f  R}  =  r  f  c!  (x)  [g*e  *e  ]  (x,  t)  dV 
t  2  J  ljkl  ~  n  kl 


V 


ij 


V 

+  2 


p(x)  p!.(x)  [g+Tfe  ]  (x,  t)  dV 
~  ij  ij  ~ 


r  T° D(2) 
l  Ce(2> 


[g^n^n]  (x,  t)  dv 


i 

+  2 


p(x)  [u.*u  ]  (x,  t)  dV 

V 

l  r  kii(~} 

-j  — £ - [g*g'*,?.*t?  ]  (x,  t)  dV 

V  o  J 


V 


[g*T  ,*e  ]  (x,  t)  dV 

ij  ij  ~ 


p(x)  [g*T)*9j  (x,  t)  dv 


V 


V 


[g*g'*q.*t>  ]  (x,  t)  dV 


1  1 


-  f  Ug*T  +  f.)*u  ]  (x,t)  dV 
«  1  j  i  J  1  1 

V 


— -  [g*(g'*q  -  h)*9 ]  (x,  t)  dV 

i  i  j  l 

V  o 

[g*T.*u.]  (x,  t)  dS 


+  f  [g*(T  -  T  )*u  ]  (x,t) 
J  i  l  l  ~ 


dS 


+  J  y~  [g*g,*Q*e]  (x,  t)  ds 
+  f  ^-[g*g'*(Q-  Q)*9]  (x,t)  ds, 

%  To 


(5.2) 
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where 


and 

Then 

6Q^{  R}  =0  over  K  (0^t<oo)  (5.4) 

if  and  only  if  R  is  a  solution  of  the  mixed  problem 

Proof.  Let  R  =  fu,  je,  T,  "9*,  T),  c[,  *?]  e  K,  from  which  it 
follows  that  R  +  X  R  e  K  for  every  scalar  X.  Then,  by  (5.  2),  (2.  7), 
(3.4),  (5.  Id),  (5.  le),  property  (b)  of  admissible  states,  the  proper¬ 
ties  of  the  convolution,  and  the  divergence  theorem, 


P!.(x)  =  —(  v  r  7  v  P.  .(x) 
L  P(x)  C  (x)  ij - 


(5.  3a) 


C'  ,(x)  =  c.  (x)  +  3.  (x)  P>  (x) 
ijkl  -  ljkl  -  ij  -  kl- 


CS.  3b) 


6RntfRl  =  J  c;jkl  -ki  -  p  v  -  v*V  <2.  t>  dv 

T  X] 

+  |  p(x)  [g*(-^4-  -  ej.  e.  -  e ) *tt]  (x,  t)  dv 

v  e  1J  1J 

-  [  [  g*  g 1  * (k .  t?  +  q.)*t?  ]  (x,  t)  dV 

J  1  it  i  i  i 

V  o 

-  [  [(g*T..  .  +  f.  -  PU  )*u.]  (x,  t)  dv 

ij.  J  i  ii- 

-  f  TfL  [g*{g'*q.  .  -  h  +  p  T  Tl)*e]  (x,t)  dV 

J  i  1,1  o  ~ 


V  o 


V 


[g*(u  -  e..)*T. .]  (x,  t)  dV 

(i.j)  ij  iJ  ~ 


(5.  5) 


+  ^  -=~  [g*g'*(8, ,  -^.)*q.]  (x,  t)  dV 

J  X  111 


v  ° 


+  [  [gs:e(u.  -  u.)*T.]  (x,  t)  dS 

J  111 
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+ 


[g*(T. 


t  (x,  t)  ds 
1  1  ~- 


+ 


+ 


J-[g*g'*(0  .  0)*Q]  (x,  t)  dS 
o 

~  [g*g'*(Q  -  Q)*9]  (x,t)  dS 

O 


(0  S  t  <  »). 


First  suppose  R  is  a  solution  to  the  mixed  problem.  Then,  by  virtue 
of  Theorem  4.  3,  (5.  5)  yields 


R}  =  0  (  0  =  t  <  ®>)  fo  re  very  ReK,  (5.6) 

which  implies  (5.  4) 

Now  examine  the  ’’only  if11  portion.  It  must  be  shown  that 
ReK  is  a  solution  to  the  mixed  problem  whenever  (5.  6)  holds. 
Choose  R  =  [ul,  0,  0,  0,  0,  0,  0],  and  let  u!,  together  with  all  of 
its  space  derivatives,  vanish  on  S  x  [0,°°).  Then,  from  (5.  1),  (5.  5), 
(5.  6),  follows 


J  [(g*T.  .  +f  -  PU  )*u]  (x,  t)  dV  =  0  (OS  t<  oo),  (5  7) 

V  iJ.J  i  l  l  ~ 

_  1,2 

and  (5.  7)  must  hold  for  every  u  e  C  with  the  foregoing  proper¬ 
ties.  But  this  fact,  together  with  Lemma  2.  1  implies  the  validity  of 
(4.4).  Next  let  R  =  [u,  0,  0,  0;  0,  0,  0],  but  this  time  require  only 
that  u  vanish  on  S  x[0,®).  Then  (5.  5),  (5.6),  (4.4),  and  Lemma 
2.  2  yield 

g*(T.  -  T.)  =  0  on  S  x  [0fo)f  (5.  8) 

1  1  Cj 

and  this  result,  because  of  (2.5),  (5.  la),  implies  (3,8b).  Now  let 
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R  =  [  0,  0,  0,  0  ,  0,  0,  0] ,  and  let  0,  together  with  all  of  its  space 
derivatives,  vanish  on  Sx[0,®).  Then,  from  (5.1),  (5.5),  (5.6) 
follows 


Y"[g*(g,*qi  -  h  +  D  r|)*0]  (x,  t)  dV  =  0 

V  o 


(5.9) 


(0St<«)f 


and  (5.  9)  must  hold  for  every  0  e  C  with  the  aforementioned 
properties.  This,  together  with  Lemma  2.  1,  implies  (4.  5)  Then 
let  R  =  [  0,  0,  0,  0 ,  0,  0,  0] ,  but  require  only  that  0  vanish  on 
$  x  [  0,  oo ).  Then  (5.5),  (5.6),  (4.  5),  and  Lemma  2.  2  imply 


^~g*g'*(Q  -  Q)  =  0  on  ^  x  [<>,•),  (5  10) 

O 

and  this  result,  with  (2  5),  (5.  la),  yields  (3.  8d),  Now  let 
R  =  [0,  0,  0,  0,rf,  0,  0],  and  suppose  rf  and  all  its  space  derivatives 
vanish  on  S  x  [0,«)  By  (5.  5),  (5.  6),  and  Lemma  2.  1 


V 

c 


p.'.  e. 
ij  ij 


-  9 


0  on  V  x  [  0,  ®). 


(5  11) 


This,  together  with  (2.  5),  (5,  la),  (5  3a), implies  (3.  5c).  Now  let 
R  =  [0,  e,  0,  0,  0,  0,  0],  and  suppose  ]e  and  all  of  its  space  derivatives 
vanish  on  S  x  [0,  «).  Since  R  and  R  are  admissible,  T.  and  *e 
are  symmetric.  Thus  (3.  4a),  (3.  4f),  (5.  5),  (5.  6),  and  Lemma  2.  1 
yield 


p!  .  T1  -  T.  )  =  0 
1J  1J 


on 


Vx[0,«),  (5.12) 
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and  this  fact,  in  view  of  (2  5),  (5  la),  (3  5c),  (5  3),  implies  (3,  5a). 
Then  let  R  =  [0,  0,  0,  0,  0,  0 ,j?] ,  and  suppose  1)  and  all  of  its  space 
derivatives  vanish  on  S  x  [  0,  oo).  Using  (5.  5),  (5.  6),  and  Lemma  2.  1 
yields 

g*g'*(k.j  +  q.)  =  0  on  V  x  [  0,  =°),  (5.  13) 

and  this  result,  with  (2.  5),  (5.  la),  implies  (3,  5b).  Next  select 
R  =  [0,  0  ,T,  0,  0,  0,  O],  where  T  and  all  of  its  space  derivatives 
vanish  on  S  x  [0,o°),  and  using  (5  5),  (5  6),  the  symmetry  of  T  and 
]e,  Lemma  2.  1,  (2.  5),  and  (5.  la),  verify  that  (3.  la)  holds  Then  let 
H  -  [0,  0,T,  0,  0,  0,  0],  but  now  assume  only  that  T  and  its  space 
derivatives  vanish  on  x  [0,oo)#  Then  (5.  5),  (5.  6),  (5  Id),  (3.  la), 
Lemma  2.  3,  (5.  la),  and  (2.  5)  imply  (3.  8a)  Then  let 
H  =  [0,  0,  0,  0,  0  0],  where  2  anc^  aH  its  space  derivatives  are 

assumed  to  vanish  on  S  x  [O,00)  Use  (5.  5),  (5.  6),  Lemma  2.  1,  (2.  5), 
and  (5.  la)  to  show  that  (3.  lb)  holds.  Finally,  let  R  =  [0,  0,  0,  0,  0,  c[,  0], 
but  insist  only  that  2  and  its  space  derivatives  vanish  on  x  [0,<»). 
Then  (5.  5),  (5.  6),  (5.  le),  (3.  lb),  Corollary  2.  4,  (5.  la),  and  (2.  5) 
imply  (3.  8c),  Therefore  R  satisfies  (3.  1),  (3.  5),  (4  4),  (4.  5), 

(3.  8),  and  hence  from  Theorem  4.  2  R  is  seen  to  be  a  solution  of 
the  mixed  problem,  which  completes  the  proof, 

J3LL  Tjye^mye^m_  Let  K  be  the  set  of  all  admissible  states 
that  meet  the  strain-displacement  and  thermal  gradient-temperature 
relations  (3.  1).  Let  R  =  [u,  ,  0  ,  T),  C[,$]  €  K  and  for  each  te  [0,<») 

define  the  functional  { *  }  on  K  by 

©Jr}  =  J  [g*T.  *e..]  (x,t)  dV 
t  Jv  Ij  Ij  - 

-  J  P(*)  [g*r|*9]  (x,  t)  dV 
V 
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where  ccj. 

ij 


i(W~,(g*TiAi1(5’t)  dv 


-  J  p(x)  a!.(x)  [g*T|*T..]  (x,  t)  dV 


V 


ij 


ij 


l  T  p(x) 

+  J  J  '  c~~'i  (x.*)  dV 

+  "pJ  P(x)  tu  *u  J  (x»  t)  dV 

2V  ii~ 


+  -jJ  [ g':< g ' q . 'Sc^  J  (x,  t)  dv 

Vo  11 


1  (2) 

+  —  J  -  [g^g'^q.^qj  (X,  t)  dv 


2  J  T 
V  o 


(5  14) 


-  J  [f.*u  ]  (x,  t)  dV  +  J  -r-[g*h*9]  (x,  t)  dV 
V  1  1  V  o 


+  J  [g*(u.  -  u.)*T.  ]  (x,  t)  dS 

q  111^* 

bl 


-  J  [g*T.*u.]  (x,t)  dS 


+  J  7“[g*g'*(0  -  0)*Q]  (x,  t)  dS 
Jl  ° 


-  J  ~[g*g'*Q*0]  (x,t)  dS, 
^T° 


(x)  and  X!  ,,(x)  are  given  by  (4.9).  Then 
—  ijkl  ~ 


60^_{R}  =  0  over  K  (0  =  t  <  ®>) 


(5.  15) 


if  and  only  if  R  is  a  solution  of  the  mixed  problem. 
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Proof,  Let  R  =  [u,  ]e,  T  ."ST,  rf ,  c[;  7J  be  an  admissible  state, 
and  suppose  in  addition  that  R  +  \  H  e  K  for  every  scalar  This 

latter  condition  is  equivalent  to  the  requirement  that  R  e  K.  Then, 
from  (5.  14),  (2,  7),  (3,  1),  (3.  4),  the  properties  of  the  convolution, 
and  the  divergence  theorem,  follows 


6^0t{R] 


J  [g*(< 


V 


ij 


~  y  1  t 

ijkl  kl 


pa!,  t]  )*  t.  ]  (x,  t)  dV 
ij  iJ  ~ 


+ 


J  p(x)  [g*( 

V 


V 

C 


T  . 

IJ 


dV 


+  J  +  X  q.)*q.]  (x.t)  dV 

V  o  1  1J  J  1 

f  [(g*T.  +  f  -  p  u.)#ll  ]  (x,  t)  dV 

ir  1J.J  1  1  1  - 


-  f  7f-[g*(g'*q.  -  h  +p  T  T| (x,  t)  dV  (5.  16) 

VTo 

+  J  [g*{u  -  u.)*T  ]  (x,  t)  dS 
S1 

+  J  [g*(T  -  T.)*n.]  (x,t)  dS 
S2  111 

+  J  4-[g*g'*(e  -  0)*q]  (x,t)  ds 

° 

+  r  ~  [g*g'*(Q  -  Q)*e]  (x,  t)  ds,  o  ^  t  <  oo. 

To 


If  R  is  a  solution  of  the  mixed  problem,  then  (5  16),  because  of 
Theorem  4.  3,  yields 


29 


6^0t(R}  =0,  0^  t<  co,  for  every  R  e  K,  (5.  17) 


which  implies  (5.  15).  On  the  other  hand  (5.  16),  (5  17),  Lemma  2.  1, 
Lemma  2.  2,  Lemma  2.  3,  Corollary  2.  4,  (2.  5),  (3.  1),  and  Theorem 
4.  3  imply  that  R  is  a  solution  to  the  mixed  problem,  and  the  theorem 
is  proved. 

Another  form  of  vanational  principle  can  be  developed  by 
requiring  more  stringent  conditions  on  the  admissible  states. 
Analogous  to  the  notion  of  a  kinematically  admissible  state  as  used 
by  Gurtin  [5],  a  kinematically  and  thermally  admissible  state  is 
introduced  as  a  state  that  meets  the  strain-displacement  and  thermal 
gradient-temperature  relations  (3.  1),  the  constitutive  relations  and 
equation  of  state  (3.  5),  and  the  displacement  and  temperature  bound¬ 
ary  conditions  (3.  8a),  (3.  8c).  The  variational  principle  based  upon 
such  states  represents  an  extension  of  the  theorem  of  minimum 
potential  energy. 

5.J3 .  Theorem.  Let  K  be  the  set  of  all  kinematically  and 
thermally  admissible  states.  Let  R  =  [u,  e_, JT,  0 ,  T) e  K,  and  for 
each  te  [0,  qo)  define  the  functional  $  {'}  on  K  by 


§  {R}  =  yj  [g*T*e  ]  (x,  t)  dV 

t  2  y  ij  ij  - 

-  tJ  o(x)[g*n*e]  (x.t)  dV 

v 

+  \\  p(x)  [u.*u.]  (x,t)  dV 
2V  1  1  ~ 


+  L  T7-  [g*g'*q.*d  ]  (x,  t)  dV 
V  o  11 


(5.  18) 
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Then 


-  J  [f.*u.]  (x,t)  dV 

X,  1  1  ^ 


+  J  =r-  [g*h*9]  (x,  t)  dV 
V  o 

-  J  [g*T.*u.]  (x,  t)  dS 
S2 


-  J  ~[g*g'*6*9]  (x,t)  dS 

JZ 


6$  f R}  =  0  over  K,  05  t<» 


if  and  only  if  R  is  a  solution  to  the  mixed  problem. 

Proof .  Let  R  =  [  u,  je,  T  ,  "ST,  r\  ,  cpz?]  be  an  admissible 
suppose  that 


R  +  X  R  e  K  for  every  scalar  X. 


This  latter  condition  is  equivalent  to  the  requirement  that 
(3.  1),  (3.  5),  with 


and 


u.  =  0 
i 

?  =  0 


on  S  x  [  0, 00 ) 

on  J 1  x  [0,®). 


(5.  19) 

state  and 

(5.  20) 

meets 

(5. 21a) 

(5. 21b) 


Next  use  (5.  18),  (2.  7),  (3.  1),  (3.  4),  (3.  5),  (5.  21),  the  properties  of 
the  convolution,  and  the  divergence  theorem  to  verify  that 
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+  f 
J  i 


-  pu.  )*vl.]  (x,  t)  dV 
1  l  ~ 


-  h  +  p  T  T\  )*?]  (x,  t)  dY 


o 


+  J  [g*(T.  -  T.)*H.]  (x,t)  dS 

1  i  ii  — 

S2 


(5.22) 


0  —  t  <  «, 


R  is  a  solution  to  the  mixed  problem.  Then  (5.  22),  by  virtue  of 
Theorem  4.  3,  implies  (5.  19).  On  the  other  hand  (5.  19),  (5.  22), 
Lemma  2.  1,  Lemma  2.  2,  (2.  5),  (3.  1),  (3.  5),  and  Theorem  4.  3 
imply  that  R  is  a  solution  of  the  mixed  problem.  This  completes 
the  proof. 
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6.  a  riati  on  a  1  yh ajract  e  r  izati  on  of  d  rs  pi acem  e  n  ts  and 

tempe  rature . 


meant  an  ordered  pair  [u,  8]  such  that  the  vector-valued  function 
u  e  C^'  ^  and  the  scalar-valued  function  0  e  C1’  The  linear 
space  of  interest  here  will  be  the  set  of  all  admissible  displacement 
and  temperature  fields.  By  a  k in e mat i c a lly  and^thermany  a  dmi s s lb  1 e 
dis j)l  a  c  e_m  ent^a  n  d  jtem  per  a  ture  f  pe  1  d  j  s  meant  an  admissible  displace¬ 
ment  and  temperature  field  satisfying  the  boundary  conditions  (3.  8a) 
and  (3.  8c)  on  the  displacement  and  temperature.  It  follows  that,  if 
[u,  je ,  T^,  0  ,  T] , is  a  kinematically  and  thermally  admissible  state, 
then  [u,  0]  is  a  kinematically  and  thermally  admissible  displacement 
and  temperature  field.  Conversely,  the  latter  assertion  implies  the 
former  when  and  q  are  defined  through  (3  1)  and  (3.5) 

Thus  the  following  corollary  to  Theorem  5.  3  may  be  stated. 

6.  1.  Corollary.  Let  K  be  the  set  of  all  kinematically  and 
thermally  admissible  displacement  and  temperature  fields  Let 
[u,0]  e  K  and  for  each  t  e  [0;.  »)  define  the  functional  q^^_{  *  }  on  K 
through 


cpt{u,9}  = 


+ 


2  l  dV 


T  J*  p(x)  [u*u.]  (x,  t)  dV 
2V  ^  1  1  ~ 


-  j  P..(x)[g*u.  .*0](x,t)dV  (6.1) 

v  1J  1,J 

p(x)  C  (x) 

-  jJ  - — [g*0*e]  (x,t)  dV 

Y  o 
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,  k..(x) 

vj  — [g*g'*e,.*e, .]  (x, t)  dv 

V  o  J 


-  J  [f  *u  ]  (x,  t)  dV  +  J  7—  [g*h*9]  (3 
V  1  1  Vo 


-  J 


[g*T.*u.]  (x,  t)  dS 


J  Y”[g*g'*Q*e]  (x,t)  as. 
A 2 


Then 


6cp^{u,  0}  =  0  over  K  (0=  t<  <») 


,  t)  dV 


(6  2) 


if  and  only  if  u  and  0  are  displacement  and  temperature  fields 
corresponding  to  a  solution  of  the  mixed  problem  of  coupled  thermo¬ 
elasticity. 
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7.  Variational  characterization  of  stress  and  heat  flax, 


By  an  admis sible  stress  and ^heat_ flu^ jield^  is  meant  an 

ordered  pair  (/Ljcj]  such  that  the  symmetric  second-order  tensor- 

""2,0  .2,0 

valued  function  e  C  and  the  vector-valued  function^  e  C 

The  linear  space  basic  to  the  development  here  is  the  set  of  all 
admissible  stress  and  heat  flux  fields.  By  a  dynamically  and 
energetically  admissible  stress  and  heat  flux  field  is  meant  an 
admissible  stress  and  heat  flux  field  satisfying  the  boundary  condi¬ 
tions  (3.  8b)  and  (3.  8d)  on  the  traction  vector  and  normal  heat  flux 
7^1.  Theorem.  Let  K  be  the  set  of  all  dynamically  and 
energetically  admissible  stress  and  heat  flux  fields  Let 
R  =  [ T ,  q]  e  K  and  for  each  t  e  [0,®)  define  the  functional  A  f’} 
on  K  through 

A{R]  -  7J  x'  ,(x)  [t  -*T  ]  (x,t)  dV 

t  2  Jv  ljkl  ~  1J  kl  ~ 


a|.(x) 

+  J  4J -  [T.  *(h  -  g'*q  )]  (x,t)  dV 


T 

V  o 


ij 


m,  m 


.  X  (*> 

-  —  J  -  [g'*q.*q.]  (x»  t)  dV 

Z  -iT  1  _  1  J 


V  o 


+  tJ  7 7~r[g*T  -  *T  ]  (x,  t)  dV 

2  Jv  p(x)  ij,  J  im,  m  ~ 


(7  1] 


-•  —  f  - - - r - [g'*g,,;<q  *q  1  (x,t)  dV 

2  y  P(x)  cTl~'  T0  i  m,mJ  ~ 


+  f  "7-7  [f  *T..  ]  (x,  t)  dV 

Jv  P(x/  1  iJ.J  ~ 
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+  J 


p(x)C  (x)T 
V  —  T  —  O 


[g'*h*q.  .]  (x,  t)  dV 
1  »  1 


J  [u.*T.]  (x,t)  as 


r  Tr-[g'*9*Q]  (x,t)  dS, 
Jl  ° 


where  a!j(x)  and  are  given  by  (4  9).  Then 


6A^{R}  =  0  over  K  (0^t<  00 ) 


(7  2) 


if  and  only  if  JT  and  are  stress  and  heat  flux  fields  corresponding 
to  a  solution  of  the  mixed  problem  of  coupled  thermoelasticity. 

Proof.  Let  R  e  K  so  that 


T. 

=  "t. .  n.  =  0 

on 

s9 

X  [0,  °°) 

(V. 

3a) 

i 

ij  J 

2 

and 

Q  = 

-  "q  n  =  0 

on 

A 

X  [o,  °°) 

(7. 

3b) 

It  follows  that  R  +  \  13.  e  K  for  every  scalar  X.  Then,  by  (7.  1), 

(3.  4b),  (3.  4d),  (3.  4e),  (2.  7),  (7.  3),  the  properties  of  the  convolution, 
the  symmetry  of  T,  and  the  divergence  theorem, 


V,lRl 


J  {[X!jklTkl 


a! . 


(h  - 


g'*q_ 


m, 


) 

m 


-  (> 


Tr  +  t 

(im,  m  ( 


i,),J)K}(~'t|dv 
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+  J  7-g'*{[(akrkl-^-[h-«'*Vm1),i 


v  ° 


-  X..  qj-qj(x,t)  dV 

+  J  {CF<**Vjwi,'4i]*Ti}(5't,dS 


(7.4) 


+  f  t"  g  ( r^T7- <h’g  *q-  •  -a!-T-n 

V  L  ^PCT  1, 1  ij  ij 

Jl  °  T 

-  e]  *q}>  (x, t)  as,  o  st<  «, 


for  every  R  e  K.  If  R  corresponds  to  a  solution  of  the  mixed 
problem  then  (7.  4)  and  Theorem  4.  5  yield  (7  2).  On  the  other  hand 
(7.  2),  (7.  3),  (7.  4),  (3.  4),  the  symmetry  of  JT,  Lemma  2.  1,  Lemma 
2 .  3,  Corollary  2.  4,  (2.  5),  (5.  1),  and  Theorem  4  5  imply  that  R  is 
a  solution  to  the  mixed  problem.  This  completes  the  proof. 

In  closing  we  note  that  extensions  of  these  theorems  to  a 
theory  of  time -invariant  thermoviscoelasticity,  in  which  the  mech¬ 
anical  properties  are  independent  of  temperature,  are  easily  obtained. 
Because  of  the  strong  thermal  dependence  of  the  relaxation  moduli  of 
a  real  viscoelastic  material,  such  a  theory  is  felt  to  be  of  little 
significance. 
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